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edge  of  alkali  halides  and  alkaline  earth  fluorides.  This  equation  is  used  to 
calculate  the  intrinsic  absorption  coefficient  as  a  function  of  both  wavelength 
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of  technical  report  which  includes  discussions  on  the  absorption  of  lattice 
vibration,  free  carriers,  impurities  and  absorption  edges.  Temperature  and 
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1.  ABSTRACT 

Under  this  contract  a  study  was  carried  out  of  the  optical  properties  of 
selected  optical  materials  through  data  compilation,  critical  data  evaluation 
and  generation  of  recommended  values  based  on  available  data  from  worldwide 
literature.  The  selected  optical  materials  include  alkali  halides,  alkaline 
earth  halides,  silicon  and  germanium,  and  III-V  compounds.  A  major  contribu¬ 
tion  of  this  work  has  been  the  exhaustive  compilation  of  absorption  coefficient 
and  the  formulation  of  an  Urbach  type  equation  for  the  multiphonon  region  neur 
the  infrared  absorption  edge  of  alkali  halides  and  alkaline  earth  fluorides. 
This  equation  is  used  to  calculate  the  intrinsic  absorption  coefficient  as  a 
function  of  both  wavelength  and  temperature  in  the  temperature  region  of 
practical  interest  with  the  calculated  result  in  concordance  with  the  best 
contemporary  experimental  measurements  for  the  pure  samples. 

Optical  constants  of  silicon  and  germanium  were  studied  in  a  similar 
manner  as  for  the  above  mentioned  materials.  Refractive  index  data  for  silicon 
and  germanium  were  exhaustively  compiled  and  analyzed.  Recommended  values  of 
refractive  index  for  the  transparent  wavelength  region  were  generated  in  the 
range  1.2  to  14  pm  and  100-750  K  for  silicon,  and  1.9  to  16  pm  and  100-550  K 
for  germanium.  Generation  of  these  values  was  based  on  a  dispersion  equation 
which  best  fits  the  selected  data  sets  covering  wide  temperature  and  wavelength 
ranges.  Both  the  refractive  index  and  its  temperature  derivative  are  in  con¬ 
cordance  with  the  available  data.  Data  compilation  for  the  absorption  coeffi¬ 
cient  of  silicon  and  germanium  is  completed  and  is  currently  in  the  stage  of 
data  analysis  and  preparation  of  technical  report  which  inoludes  discussions  on 
the  absorption  of  lattice  vibration,  free  carriers,  impurities  and  absorption 
edges.  Temperature  and  pressure  dependences  of  energy  gap  are  also  presented. 

The  work  on  the  refractive  index  of  III-V  compounds  is  currently  at  the 
stage  of  data  extraction. 

2.  PERSONNEL  INVOLVED  WITH  RESEARCH  PROGRAM 

During  the  progress  of  the  work,  the  staff  members  who  have  contributed  to 
the  various  tasks  of  this  program  are  as  follows: 
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Dr.  H.  H.  Li 

Dr.  T.  C.  Chi 

Dr.  Y.  S.  Touloukian 

Mrs.  A.  F.  Furdyna 

Mr.  B.  Runyan 

Mr.  D.  J.  Huff 

In  addition  to  the  above  effort,  invaluable  contribution  was  made  by  the 
professional  personnel  of  the  Division  of  Scientific  Documentation  at  CINDAS  in 
the  literature  search,  data  source  identification  and  document  procurement. 

3 .  INTERACTIONS 

The  director  of  CINDAS,  Dr.  Y.  S.  Touloukian,  frequently  participated  in 
meetings  and  conferences  on  the  laser  technology.  He  has  discussed  the  recent 
development  of  high-power  laser  technology  with  the  attendees  and  gained  infor¬ 
mation  useful  to  this  research  project.  Dr.  Touloukian  has  been  constantly 
communicating  with  Dr.  A.  H.  Guenther  of  the  Air  Force  Weapons  Laboratory  at 
Eirtiand,  New  Mexico  regarding  the  subject  of  optical  constants  or  laser  window 
materials. 

In  addition  to  attending  meetings  on  the  related  subjects,  the  principal 
investigator.  Dr.  H.  H.  Li,  frequently  discusses  and  exchanges  information  with 
tne  faculty  of  the  Department  of  Physios  at  Purdue  University  on  the  optical 
and  other  related  properties  of  solids.  Results  derived  from  these  discussions 
are  helpful  in  locating  information  and  solving  problems, 

4.  NEW  DISCOVERIES,  PATENTS  OR  INVENTIONS 

There  were  no  new  discoveries,  patents  or  inventions. 

5.  STATUS  OF  RESEARCH  EFFORT 

The  objectives  of  this  research  project  is  to  provide  evaluated  data  and 
information  on  optical  properties  at  high  temperatures  of  selected  solid 
materials  of  potential  importance  in  laser  window  and  laser-hardened  material 
technologies. 


3 


The  goals  of  the  present  work  are: 

(1)  to  review  the  present  state  of  knowledge  on  the  refractive  index 
and  absorption  coefficient  of  selected  groups  of  pare  solid  sub¬ 
stances  related  to  the  Air  Force  environment  resistant  material 
technology, 

(ii)  to  compile  the  available  data  and  information, 

(iii)  to  critically  evaluate  and  analyze  the  existing  data,  and  to 
generate  recommended  values. 

Selection  of  substances  was  based  on  the  needs  of  the  Air  Force  material 
technology  and  the  recommendation  of  the  National  Material  Advisory  Board  on 
high-power  infrared-laser  windows.  The  research  program  planned  to  continue  to 
investigate  the  following  groups  of  substances: 

(i)  alkali  halides 

(ii)  alkaline  earth  halides 

(iii)  diamond-structures  materials  (Ge  and  Si) 

(iv)  III -V  compounds  (GaP,  GaAs,  GaSb,  InP,  InAs,  and  InSb) 

(v)  (ZnS,  ZnSe,  ZnTe,  CdS,  CdSe,  and  CdTe) 

(vi)  silver  halides  (AgCl,  AgBr,  and  Agl) 

(vii)  other  IR  materials  to  be  selected. 

Ve  have  completed  the  works  for  material  groups  (i)  and  (ii).  At  the  time  of 
writing  this  final  technical  report,  the  work  on  the  material  group  (iii)  is  in 
the  stage  of  data  analysis  and  technical  report  preparation.  The  work  on  the 
material  group  (iv)  is  in  the  stage  of  data  extraction. 

5.1.  Studies  on  the  Absorption  Coefficient  of  Alkali  Halides 

A  typical  absorption  spectrum  of  alkali  halide  is  shown  in  Figure  1.  The 
fundamental  transparent  region  cf  alkali  halides  is  defined  between  the  Drbach 
tail  and  multiphonon  absorption.  In  the  high  transparent  region,  the  trans¬ 
parency  of  the  material  is  limited  by  many  factors,  notably  the  crystal 
vacancies,  dislocations,  impurities,  surface  contamination,  etc.  Extrinsic 
absorption  due  to  these  origins  can  be  reduced  through  appropriate  crystal 
growing,  purification,  annealing,  and  surface  cleaning  processes. 

Investigation  of  the  Urbach  taxi  .region  may  serve  as  an  indicator  to  show 
the  extent  of  impurity  and/or  defoct  contnnis.  As  generally  observed,  the 
purer  the  sample  is,  the  more  the  validity  of  Urbach  rule  is  extended  to  the 
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Figure  1.  Absorption  Spectrum  of  a  Typical  Alkali-Halide  Crystal 
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transparent  region.  By  comparing  the  absorption  spectrum  at  the  Urbach  tail, 
the  purer  sample  is  self-revealed. 

Absorption  in  the  multiphonon  region  has  been  of  current  interest  because 
of  its  application  in  windows  for  high-power  infrared  lasers.  At  the  frequen¬ 
cies  in  this  region,  the  total  absorption  can  be  attributed  to  processes  in¬ 
volving  several  phonons,  defect  modes  due  to  impurities,  vacancies,  or  surface 
contaminations.,  A  number  of  observations  have  been  made  to  investigate  the 
frequency  and  temperature  dependence  of  the  intrinsic  multiphonon  absorption. 
It  has  been  found  that  the  exponential  dependence  of  the  absorption  coefficient 
on  frequency  holds  for  LiF,  NaCl,  KC1 ,  and  BBr  at  room  temperature,  i.e., 

-\V  v 

a  *=  a  e  0  (1) 

o 

Ve  have  found  that  this  exponential  relation  is  also  applicable  in  the  cases  of 
NaF  and  KI. 

This  exponential  relation  attracted  considerable  attention  in  the  theoret¬ 
ical  interpretation  of  such  behavior.  Two  fundamentally  different  approaches 
have  been  taken  to  calculate  the  anharmonically  induced  absorption  in  ionic 
crystals.  The  most  common  approach  requires  the  solution  of  the  harmonic- 
lattico  problem  followed  by  a  perturbation  treatment  of  the  enharmonic i ty . 
Sham  and  Sparks  and  McGill  et  al.  employ  diagrammatic  techniques  for  evaluating 
the  Green's  functions,  while  Bendow  et  al.  use  the  equation-of-motion  method. 
Calculations  using  these  methods  become  increasingly  complex  for  higher-order 
processes  and  consequent  simplification  and  approximation  are  made. 

The  results  indeed  predict  the  exponential  dependence  of  absorption  coef¬ 
ficient  on  the  frequency.  With  regard  to  the  temperature  dependence,  however, 
the  results  are  not  satisfactory.  Discrepancies  between  experimental  and 
theoretical  predictions  are  of  several  orders  of  magnitude.  In  the  case  of 
NaF,  NaCl  and  KC1,  for  example,  theory  predicts  that  temperature  dependence  of 
absorption  coefficient  follows  the  relation; 

a  -  Tn“1  (2) 

where  the  values  of  n  are  4,  6  and  7  for  NaF,  NaCl  and  KC1,  respectively,  but 
the  corresponding  experimental  values  are  2.6,  3.3  and  2.8. 


In  our  studies,  we  have  compiled  available  data  for  various  spectral  and 
temperature  regions  that  have  been  experimentally  investigated.  For  the  data 
in  the  multiphonon  region  we  found  that  the  exponential  relation  also  holds 
well  at  various  temperatures  and  that  the  curves  tend  to  converge  to  a  point  in 
t/  low  frequency  region.  This  situation  is  analogous  to  the  Drbach  rule  in 
the  ultraviolet  absorption  edge.  Although  the  experimental  temperature  depen¬ 
dence  of  absorption  coefficient  does  not  agree  with  theoretical  prediction  as 
mentioned  abr-re,  the  observed  trend  does  follow  a  power  relation:  a  *■  Ta 
except  at  low  temperatures.  Based  on  these  facts,  it  appeared  that  an  equation 
of  ' e  fo^w: 

a  ■=  a  e-*<b+V)(c-log  T>  (3) 

o 

cor  -?  \  .  ,  ■  mA  to  represent  the  frequency  and  temperature  dependence  of  absorp- 
ti  i  -oe.'f icxenl  at  temperatures  higher  than  200  K.  Indeed  we  have  found  that 
t*  w  „u.  •ei.-ca?  tl*~a  fit  the  expression  satisfactorily.  The  best-fit  values  of 
tne  adjustable  constants,  aQ,  a,  b,  and  c  are  given  in  Table  1. 

Table  1.  The  Parameters  of  Alkali  Halides 


Crystal 

-1 

a  .cm 
o 

a,  cm 

b,cm  ^ 

c 

LiF 

104.5116 

44.08 

0.002237 

5.434 

NaF 

io5-4453 

81.87 

0.00437 

5.434 

NaCl 

104*9778 

94.478 

0.00591 

5.434 

KC1 

106.4187 

230.241 

0.00700 

5.434 

KEr 

io4*  943,3 

90.245 

0.00886 

5.434 

KI 

.3105 

207.251 

0.0098 

5.434 

Analogous  to  the  Urbach  rule,  the  pair  of  constants,  ao  and  b,  defines  the 
’cross  point’  where  the  curves  of  a  versus  V  converge.  The  physical  meaning  of 
this  point  remains  to  be  interpreted.  Figures  2  and  3  provide  a  visual  com¬ 
parison  of  the  predicted  values  from  eq.  (3)  and  experimentally  determined  data 
of  sodium  chloride.  Figure  4  is  a  plot  of  the  corresponding  calculated  values 
for  a  wider  wavenumber  region  including  the  crossover  point  (b,a  ). 
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5.2.  Studies  on  the  Absorption  Coefficient  of  Alfcslins  Earth  Fluorides 

In  this  study,  ve  compiled  all  the  available  experimental  data  of  absorp¬ 
tion  ooeffieient  and  related  information  useful  for  the  Air  Force  material 
technology  and  laser  research.  In  the  data  analysis  on  the  available  data  of 
laser  interest,  ve  found  that  eq.  (3)  can  be  equally  used  to  describe  the 
absorption  coefficient  as  a  function  of  both  frequency  and  temperature  vith 

appropriately  determined  constants,  a  ,  a,  b,  and  c  for  eaoh  material.  These 

o 

constants,  listed  in  Table  2,  were  determined  based  on  data  evaluation,  corre¬ 
lation, and  synthesis.  In  the  example  of  strontium  fluoride.  Figures  5  and  ( 

A 

provide  a  visual  comparison  of  the  predicted  values  of  eq.  (3)  and  experimen¬ 
tally  determined  data.  Figure  7  is  the  corresponding  calculated  value  for  a 

wider  wavenumber  region  including  the  crossover  point  (b,a  ). 

o 

Table  2.  The  Parameters  of  Alkaline  Earth  Fluorides 


Crystal 

log  a 

a, cm 

b,cm  * 

c 

KgF2 

3.1877 

0.004184 

-377.29 

4.2901 

c.f2 

3.0834 

0.005251 

-383.92 

4.4498 

SrF2 

2.91(9 

0.00(577 

-3(0. 5( 

4.3288 

b.f2 

2.8507 

0.007349 

-327.(7 

4.2873 

As  a  counterpart  of  the  Urbach  rule  for  the  uv  absorption  edge,  ve  have  estab¬ 
lished  the  expression,  eq.  (3),  for  the  infrared  absorption  edge  of  alkaline 
earth  fluorides.  These  expressions  are  of  the  same  form  and  the  parameters  in 
the  corresponding  equations  are  similar.  Compared  vith  Deutseh's  (eq.  1) 
expression,  ve  have  extended  the  dependence  of  the  absorption  coefficient  to 
include  the  temperature  in  addition  to  frequency.  It  is  vorthvhile  to  point 
out  that  the  equation  formulated  in  this  work  is  of  the  same  type  as  it  vas 
found  for  alkali  halides.  Correlation  techniques  used  fox  alkali  halides  are 
equally  applicable  in  the  present  vork. 

5.3.  Studies  on  the  Optical  Constants  of  Silicon  and  Germanium 

The  vork  on  the  optical  constants  of  silicon  and  germanium  is  in  the  stage 

of  data  analysis  and  technical  report  preparation.  Refractive  index  data  are 
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Absorption  Coefficient *  cm 
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•Tillable  in  the  fundamental  transparent  wavelength  region  and  ere  represented 
by  the  following  dispersion  equation  through  data  analysis: 

n2(A,T)  »  s(T)  +  (A  +  A,T  +  A.T2)  <4> 

A2  0  1  2 


where 


X  ■  wavelength  in  units  of  pm, 
T  ■  temperature  in  nnits  of  K, 


-3AL(T)/L„ 


L(T)  -  e 


»  length  of  the  sample  at  temperature  T, 


a(T)  “  statio  dieleotrio  constant  at  temperature  T,  and 
Ao,  Aj,  and  A^  are  constants  for  a  given  material. 

In  the  case  of  silicon,  the  parameters  are: 

«(T)  -  11.4445  +  2.7739  x  10"4T  +  1.7050  x  lO-6!2  -  8.1347  x  lO"10!3, 
Afl  -  0.8948,  k±  -  4.3977  x  10~4,  Aj  -  7.3835  x  10~8. 


-  -0.021  -  4.149  x  10“7T  -  4.620  x  lO^T2  +  1.482  x  10-11!3 

L293 

(for  temperature  region  20-293  K) 


-  -0.071  +  1.887  x  1P~^T  1.934  x  10 "9T*  -  4.544  x  10-13!3 

L293 


(for  tenperati 


igion  293-1600  K) 


In  the  case  of  germanium,  the  parameters  are: 


e(T)  -  15.2892  +  1.4549  x  10~3T  +  3.5078  x  lO^T1  -  1.2071  x  10'  ■V. 
Aq  -  2.5381,  Aj  -  1.8260  x  10”3,  Aj  -  2.8888  x  10~6, 


MAH  .  _j 


0.089  +  2.626  x  lO^T-lOO)  +  1.463  x  10“®(T-100)2  - 


2.221  x  10_11(T-100>3 


(for  temperatures  100  <  T  <  293), 


and 
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**  5.790  -*  *  1.768  x  10"9<T-293)2  - 

L293 

4.562  x  10-13(T-293)3  (for  temperatures  293  <  T  <  1200). 

Figure*  8  sad  9  show  the  calculated  refractive  index  and  its  temperature 
coefficient*  dn/dT*  of  silicon*  the  latter  is  calculated  from  the  first 
temperature  derivative  of  eq.  (4) . 

A  large  amount  of  data  on  the  absorption  coefficient  of  silicon  and 
germanium  are  oompiled  and  organised.  Shown  in  Figures  10  and  11  are 
collective  plots  of  all  compiled  data  for  silioon  and  germanium*  respectively. 
In  each  of  the  plots  there  are  about  5000  data  points  obtained  from  various 
data  sources.  Similarity  of  these  two  plots  evidences  that  the  absorption 
mechanisms  in  silioon  and  germanium  are  fundamentally  the  same  as  indicated  in 
the  plots.  The  absorption  can  be  classified  and  discussed  in  terms  of  three 
origins*  namely  the  electron  excitation,  the  carrier  absorption*  and  the 
lattice  absorption. 

a.  Electron  Excitation 

In  the  transparent  region,  the  absorption  is  usually  low.  As  the  photon 
frequency  is  increased  approaohing  the  high  frequency  edge*  a  threshold  point 
is  reaohed  where  absorption  is  abruptly  increased  to  several  orders  of  magni¬ 
tude  within  a  narrow  frequency  range  as  a  result  of  excitation  of  electrons 
from  valence  band  to  conduction  band  across  the  energy  gap.  Although  the 
threshold  of  continuous  optical  absorption  facilitates  a  means  for  the  deter¬ 
mination  of  energy  gap*  difficulties  are  encountered  because  the  process  of 
such  excitation  in  silioon  and  germanium  is  indirect*  the  participation  of 
phonons  in  the  indirect  process  must  be  considered.  A  considerable  amount  of 
work  was  dedicated  in  this  particular  problem  for  the  determination  of  acourate 
energy  gap  rs  a  function  of  tempercture  and  pressure. 

b.  Carrier  Absorption 

In  the  insulators,  optical  absorption  due  to  carriers  is  negligible.  Fhile 
in  the  metals,  optical  absorption  due  to  free  electrons  is  extremely  high) 
moderate  variation  of  the  free  electron  concentration  is  not  sensitively  re¬ 
flected  in  the  change  of  absorption  level.  Only  in  the  semiconductors*  optical 
absorption  is  very  much  influenced  by  the  carrier  concentration.  Since  the 
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FIGURE  8.  RECOMMENDED  n-X-T  DIAGRAM  OF  SILICON 


FIGURE  9.  RECOMMENDED  dn/dT-A-T  DIAGRAM  OF  SILICON 
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10.  Absorption  Coefficient  of  Silicon 


Free  carrier  absorption 


Figure  11.  Absorption  Coefficient  of  Germanium 
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type  and  concentration  of  carrier  can  be  introduced  in  a  controllable  manner, 
the  effects  of  carriers  on  the  optical  absorption  can  be  stndied.  By  use  of 
the  lav  of  mass  action,  one  can  introduce  donors  so  that  the  carriers  in  the 
sample  is  predominantly  electrons.  Absorption  dne  to  free  electrons  as  a  func¬ 
tion  of  wavelength,  temperature  and  other  parameters  can  be  differentiated. 
Likewise,  introduction  of  acceptors  makes  sample  of  hole  predominating?  thus 
absorption  due  to  holes  can  be  observed.  At  high  temperatures,  the  optical 
absorption  of  a  low  impurity  germanium  sample  behaves  closely  to  that  of  an 
intrinsic  sample  because  the  thermally  excited  electron-hole  pairs  predominate 
carriers  from  impurities.  For  heavily  doped  silicon  and  germanium,  the 
carriers  are  essentially  thermally  ionized  from  impurity  atoms.  There  is  a 
temperature  region  called  exhaustion  region  in  which  all  the  available  impurity 
atoms  are  thermally  ionized.  Since  the  impurity  concentration  is  far  more  than 
the  carrier  concentration  that  are  thermally  excited  from  valence  band  to 
conduction  band  of  the  host  crystal,  the  carrier  concentration  remains  as  a 
constant  in  the  exhaustion  region.  The  existence  of  such  a  region  offers  a 
practical  way  to  investigate  the  effect  on  the  absorption  of  parameters  other 
than  the  carrier  concentration. 

c .  Lattice  Vibration  Absorption 

In  a  pure  specimen,  the  radiation  due  to  lattice  vibration  is  conspicuous 
at  low  temperatures  while  that  due  to  thermally  excited  carriers  is  predominant 
at  high  temperatures.  In  a  doped  specimen,  the  absorption  due  to  carriers  is 
predominant  even  at  low  temperatures.  It  is  therefore  clear  that  absorption 
spectrum  due  to  lattice  vibration  can  only  be  observed  when  a  highly  purified 
specimen  is  used. 

There  have  been  a  number  of  recent  studies  on  how  the  infrared  absorption 
drops  off  as  the  frequency  becomes  much  greater  than  the  fundamental  lattice 
frequencies.  These  studies  have  been  mainly  carried  out  for  laser  window 
materials.  In  the  highly  purified  samples,  the  absorption  coefficient  exhibits 
an  exponential  fall-off  over  two  to  four  decades  to  the  lowest  values  of  a  that 
can  be  measured.  Whether  this  exponential  absorption  tail  is  characteristic  of 
all  other  classes  of  materials  is  not  known. 

All  the  absorption  mechanisms  will  be  discussed  in  the  technical  report 
which  is  currently  in  preparation.  The  corresponding  data  of  the  mechanisms 
will  be  analyzed  so  that  either  the  most  probable  values  or  the  best  available 
data  for  the  specified  conditions  are  recommended. 


5.4.  Studies  on  the  Refractive  Index  of  III-V  Compounds 


The  work  on  the  refractive  Index  of  III-V  compounds  is  currently  in  the 
process  of  data  extraction. 
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